Non-uniformity of a planar polarizer as a reason of spin-transfer induced vortex 

oscillations at zero field 
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We discuss a possible mechanism of the spin-transfer induced oscillations of a vortex in the 
free layer of spin-valve nanostructures, in which the polarizer layer has a planar magnetization. 
We demonstrate that if such planar polarizer is essentially non-uniform, steady gyrotropic vortex 
motion with large amplitude can be excited. The best excitation efficiency is obtained for a circular 
magnetization distribution in the polarizer. In this configuration, the conditions for the onset of the 
oscillations depend on the vortex chirality but not on the direction of its core. 
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Sub-GHz dynamics of magnetic vortices induced by the 
spin transfer effect observed recently in nanopillars and 
nanocontacts [ires] have raised a strong interest. The as- 
sociated microwave emissions in such Spin Transfer Vor- 
tex Oscillators (STVOs) occur without any external mag- 
netic field and at low current densities, together with 
large powers (in the nW range) and narrow linewidths 
(< 1 MHz) comparatively to single-domain spin trans- 
fer nanooscillators. In arrays of nanocontacts, a coher- 
ent motion of coupled vortex dynamics generated by the 
spin transfer, resulting in a significant improvement of 
the quality factor of the devices has been recently demon- 
strated [7j. This makes STVOs of considerable practical 
interest for new set of applications in microwave tech- 
nologies or magnetic memories. 

STVO consists of at least two magnetic layers sepa- 
rated by a non-magnetic spacer. One of the magnetic 
layers (the free layer) has a vortex that can be excited 
by the spin transfer, while the second magnetic layer is 
used as a spin polarizer for the current. So far, the theo- 
retical analysis of the spin transfer vortex dynamics has 
only considered the approximation of a fixed uniformly 
magnetized polarizer [3l. l8Hll|. For such polarizers, only 
the component of the spin polarization that is perpen- 
dicular to the plane can induce steady vortex precession 
[12\. However, in many recent experiments spin transfer 
driven vortex oscillations have been detected at zero or 
in-plane bias magnetic field, in nanopillar or nanocontact 
STVOs for which the magnetization of the polarizer nat- 
urally lies in the plane 0,1343 • The onset of a small per- 
pendicular component of the spin polarization due to the 
magnetization dynamics in the polarizer can be assumed, 
but such contribution can not be sufficient to account for 
the large amplitude vortex excitations. 

In our present work we consider another mechanism for 
the vortex excitation, which is specifically related to the 
STVO having planar magnetization distribution within 
the polarizer. First we present an analytical model for 
the vortex dynamics in a circular spin-valve nanopillar. 
The free magnetic layer of the spin valve is in a centered 



vortex state. The second magnetic layer (the polarizer) is 
magnetized in the layer plane, which leads to an in-plane 
spin polarization. To clarify our analysis, we disregard 
the stray magnetic field emitted by the polarizer and as- 
sume that it is fixed. The current flow is assumed to be 
uniform through the pillar diameter. The spin transfer 
torque [l3( is calculated using (aJ/M s )M x (M x p), 
where J is the current density, M s is the magnetization 
of saturation, M is the magnetization vector, p is the 
spin polarization vector, and a represents the efficiency 
of the spin transfer: a = his/(2 \e\ LM S ), v is the spin 
polarization of the current, e is the electron charge, L 
the layer thickness. 

We have recently suggested to use the energy dissipa- 
tion approach to derive the generalized Thiele equation 
for the spin current- induced vortex core motion [8(: 
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Here X is the vortex core position, the gyrovector is given 
by G = -Ge z , with G = 2ttM s L/j and W(X) is the 
potential energy of the moving vortex [l4| . The damping 
constant D is given by D = arjG, where the factor r\ is 
of the order of unity [15j. The last term F$t is the spin 
transfer force. 

We derive the terms of Eq. Q] using the standard two- 
vortices ansatz (TVA) for the in-plane magnetization 
components of the moving vortex [14j . The out-of-plane 
magnetization component M s cos 9 can be approximated 
by a bell-shaped ansatz suggested by Usov et al. fl6j |: 



2Ptan-! (p/b) (p < b) 
Ptt/2 (p > b) 
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where b is the core radius and P is the polarity of the 
vortex core (P = +1 if the polarity is along the z axis 
and P = —1 if it is opposite). We consider in these 
calculations a steady-state circular motion of the vortex 
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where a is the orbit radius, ui is the gyration frequency. 
The vortex energy W(~K) has two main contributions: 
the magnetostatic energy, which originates from the vol- 
ume magnetic charges arising from a shifted vortex [l4j . 
and the contribution of the Oersted field || [ItJ • The last 
two terms of Eq. ([T]) can be calculated using the energy 
dissipation function W — fw(r) dV, w is the energy 
density at point r. We use w = (SE/S8) 9 + (8E/6ip) <p 
Q, where 9 and ip are respectively the polar and az- 
imuth angles of the magnetization vector M. 5E/56 and 
5E/5ip are taken from the LLG equation. This gives us 
the current-dependent contribution to the energy dissi- 
pation density wst for the planar polarizer: 

wst = M s aJ (p x sin <p — p y cos <p) 9 

+ sin cos (p x cos <p + p y sin tp) <p] (4) 

where p x and p y are the x- and y-components of the 
local spin polarization p. The right-hand side of Eq. 
(|3J) vanishes outside the core, from which we find that 
for the planar polarizer the spin torque excites only the 
vortex core. The spin-transfer force is given by F$t = 
d (/ w ST dV) /<9X. Using Eq. © and ©, we find: 



F ST = TrM s LbPaJ (p (X) • e x ) e 



XI 



(5) 



in which the terms proportional to b 2 have been ne- 
glected. The damping force Fdamp = — D1L can be ob- 
tained by treating similarly the contribution to W pro- 
portional to the Gilbert damping a Q: 
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The vortex energy gain, given by the dot product 
(Fst + Fdamp) • X, should average to zero in each cycle 
of the steady core gyration. This leads to the following 
general condition for the onset of the steady vortex pre- 
cession for arbitrary magnetization distribution p (a, x) 
in a planar polarizer: 
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If the magnetization in the planar polarizer is uniform, 
i.e. p (x) = const, the spin transfer torque contributes 
positively to the energy gain for one semi-cycle of the 
vortex motion, (p ■ e x ) > 0, but it contributes negatively 
and with the same amplitude for the other semi-cycle, 
(p • e_ x ) < 0. Therefore a uniform planar polarizer 
should not excite the steady vortex motion. 

Nevertheless the vortex motion can be excited if p (x) 
is a non-trivial function. We conclude from Eq. (J7J that 
the planar polarizer is the most efficient when p has a 
circular distribution, that is (p • e x ) = ±1 for each Xi 
corresponding to a centered vortex in the polarizer layer 
@. For such a circular planar polarizer, we find from Eq. 



an expression for the radius of the vortex core orbit 
a as a function of the current density J: 
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where C v and C po i are the chiralities of the vortex respec- 
tively in the free and polarizer layer. From Eq. (|8"|). we 
conclude that, at a given current sign, the onset of the 
oscillations (a > 0), is not sensitive to the vortex core 
polarity P but depends on the chirality of the vortex 
C v . This feature is different from the conditions for spin 
transfer vortex excitations by a perpendicular polarizer 
as discussed below. 

We now compare our analytical results to numerical 
micromagnetic simulations, which have been performed 
for a nanopillar of 200 nm in diameter with a free NiFe 
layer of a thickness of 15 nm. We use the following mag- 
netic parameters: M s = 800 emu/cm 3 , A = 1.3 x 10~ 6 
erg/cm and a = 0.01 (values for NiFe) and a mesh with 
the cell size 2x2x5 nm 3 . The spin current polarization 
is taken to be v = 0.3. The micromagnetic simulations 
are performed by numerical integration of the LLG equa- 
tion using our micromagnetic code SpinPM based on the 
forth order Runge-Kutta method with an adaptive time- 
step control for the time integration. 
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FIG. 1: (color online) Numerical result for spin current- 
induced vortex gyration for the vortex polarizers, shown at 
the bottom (colorscale shows z-component of the polariza- 
tion). Graph: frequency / (left scale; circles: circular polar- 
izer, squares: vortex polarizer) and radius of the vortex core 
orbit a (right scale; down triangles: circular polarizer, up tri- 
angles: vortex polarizer) as function of the current density J. 
Solid line shows prediction for a( J) by Eq. [8] 

We assume that the chirality of the vortex in the free 
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layer C v is set by the symmetry of the Oersted field [l8[ . 
The polarizer layer has an artificially designed perfectly 
circular magnetization distribution, with |(p, e x )| = 1 for 
each point; thus it lies in-plane even in the disc center. 
We refer to such an idealized planar configuration as a 
circular polarizer. Its chirality is C po i . The positive cur- 
rent is defined as a flow of electrons from the free layer 
to the polarizer. 

For C v Cpoi = 1, we see that the vortex motion is ex- 
cited at positive currents, and for C v C po i = — 1 it is ex- 
cited at negative currents. We also find that results of 
the simulation are identical for both polarities of the vor- 
tex core. As shown in Fig. ([1]), a steady circular motion 
is observed for current densities larger than Jci =5 x 
10 6 A/cm 2 and smaller than J C i = 1.0 x 10 8 A/cm 2 . 
For Jci < J < Jc2, the frequency increases with J from 
0.66 GHz up to 0.90 GHz. The radius of the vortex gyra- 
tion also increases with the current, with the maximum 
value of about 60 nm. For currents densities larger than 
Jc2, the vortex polarity is periodically switching. At 
each switching event, the vortex core starts to move in 
the opposite direction [ljj . However since the spin cur- 
rent provides the excitation for both polarities on equal 
basis, the core is again accelerated by the spin torque un- 
til its velocity reaches the critical value required for the 
reversal [lllllj]. 

These numerical results are in very good agreement 
with the analytical conclusions. First, they confirm the 
possibility to excite a spin-transfer induced vortex motion 
for a purely planar polarizer. We find that the current 
sign for which the motion can be excited depends on the 
product C v Cp i but not on the core polarity P, in agree- 
ment to the predictions of Eq. (J5|). In Fig. Q] we plot as 
a solid line the prediction for the orbit radius calculated 
from Eq. (jHJ) which is in reasonable agreement with the 
numerical result [22j ]. Some quantitative difference be- 
tween the analytical and the numerical calculations can 
be ascribed to the considerable deviation of the magne- 
tization distribution from the TVA for a > b, as analyses 
of the magnetization distributions reveal. 

Both the analytical and numerical results are qualita- 
tively different from what has been found for the case of 
the vortex excitation by a uniform out-of-plane polariza- 
tion. In that case, the spin-transfer force originates from 
the out-of-the-core region of the vortex (in contrast to the 
case of the planar polarizer, for which the force originates 
from the core). It is given by F^ p = nM s p z LaJae x 
[1, El: where p z is the out-of-plane spin polarization. It 
does not depend on the vortex polarity P. Therefore, 
at a given current sign, it can excite the vortex motion 
(i.e. it can be oppositely directed to the P-dependent 
damping force Fd am p) for only one vortex polarity. In 
the steady oscillation regime, the radius of the vortex or- 
bit a is inversely proportional to small non-linear terms 
in Fg^ p and Fdamp @, EH- Due to this reason a can 
increase very rapidly with the current for J > Jqi as we 



found recently 8]. In contrast to it, for the planar po- 
larizer for J > Jc i, a (J) is proportional to the principle 
values of the forces , see Eq. ((8} ; this is the reason of the 
considerably slower dependence of a on J found in the 
simulations. 

The closest experimental situation to the idealized cir- 
cular polarizer is the polarizer in the vortex state [6j. 
For such polarizer, an additional contribution to the spin 
transfer term might arise from the out-of-plane core re- 
gion. We performed micromagnetic simulations for such 
configuration which we refer to as a vortex polarizer. We 
find, see FigfTJ that the oscillation frequency is practi- 
cally identical to the circular polarizer case, while the 
radius of the trajectory is only slightly shifted towards 
higher values. This result allows us to rule out the con- 
tribution of the vortex core in the planar polarizing layer 
as the major source of polarization for the spin transfer 
force. 

As a conclusion, we demonstrate both analytically and 
by micromagnetic simulations that a non-uniform planar 
polarizer can induce spin-transfer vortex oscillations. In 
the case of an ideal circular polarizer, we have derived 
the conditions for sustained vortex precession in the free 
layer as a function of the current signs and the respec- 
tive vortex chiralities. Important conclusion is that the 
condition for spin transfer precession does not depend on 
the core polarity. In addition, at large currents, a mul- 
tiple back and forth core switching during the motion is 
predicted. 
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